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ABSTRACT 

We analyze the far-ultraviolet (FUV) spectra of 33 classical T Tauri stars (CTTS), 
including 20 new spectra obtained with the Advanced Camera for Surveys Solar Blind 
Channel (ACS/SBC) on the Hubble Space Telescope. Of the sources, 28 are in the ~1 
Myr old Taurus- Auriga complex or Orion Molecular Cloud, 4 in the 8-10 Myr old Orion 
OBla complex and one, TW Hya, in the 10 Myr old TW Hydrae Association. We also 
obtained FUV ACS/SBC spectra of 10 non-accreting sources surrounded by debris disks 
with ages between 10 and 125 Myr. We use a feature in the FUV spectra due mostly 
to electron impact excitation of H2 to study the evolution of the gas in the inner disk. 
We find that the H2 feature is absent in non-accreting sources, but is detected in the 
spectra of CTTS and correlates with accretion luminosity. Since all young stars have 
active chromospheres which produce strong X-ray and UV emission capable of exciting 
H2 in the disk, the fact that the non-accreting sources show no H2 emission implies that 
the H2 gas in the inner disk has dissipated in the non-accreting sources, although dust 
(and possibly gas) remains at larger radii. Using the flux at 1600 A, we estimate that 
the column density of H2 left in the inner regions of the debris disks in our sample is 
less than ~ 3 x 10~^ g cm~^, nine orders of magnitude below the surface density of the 
minimum mass solar nebula at 1 AU. 
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quence 



1. Introduction 



Gas comprises 99% of the mass of primordial disks. As time increases, it is accreted onto the 
star, formed into planets, and lost by photoevaporation, leaving behind a debris disk, in which 
most of the mass is locked into planets and other solid bodies traced by secondary dust arising 
from collisions. Although the general outline of this process is agreed upon, many specific questions 
remain unanswered, mainly because the gas is difficult to observe. As a result, only ~1% of the disk 
mass, the dust, has been used as a probe of the disk evolution. However, although interconnected, 
the evolution of gas and dust may take different paths (jPascucci et al.ll2009l ) , making observations of 
the gas itself necessary to understand these processes. Of particular importance are observations of 
the gas in the inner disk, because it sets the chemical and physical conditions for planet formation. 
The bulk of the gas in these cold disks is in H2, which lacks a permanent dipole component, so the 
pure rotational and rovibrational lines are weak. Nonetheless, extensive surveys of these lines in 
primordial disks have been carried out (Bary et al. 2008; Bitner et al. 2008, and references therein), 
and they have been detected in a handful of objects. Searches using less abundant molecules 
have also succeeded and provided information on the gas in the inner regi on of gas-rich disks 
(jCarr NajitJbood : Isalyk et ahlbood : IPascucci et al.ll2009l : iNajita et aklbood ). Gas has also been 
searched for in disks of more evolved sources which are no longer accreting, within the age range 
when the trans i tion f rom primordial to debris is supposed to happen, ~5 - 20 Myr. In particular, 
Pascucci et al.l (|2006l ) looked for H2 in the disks of several non-accreting sources and found that 
the amount of gas still present at 5 - 20 Myr is not large enough to form the gas giant planets at 
that time. This observation agrees with results indicating that the am ount of hot gas in dis ks of 
non-accreting sources is decreased when compared to accreting sources (jCarmona et al.l 120071 ). 



UV observations are very promising for detecting the gas. The strong stellar Lya radiation 
bathes the UV thin regions of the circumstellar material and, as long as the H2 has a temperature 
of a few thousand degrees, the line excites electrons to u pper electronic state s, which produces a 
plethora of emission lines in the UV when they de-excite (jHerczeg et al.ll2006l . H06, and references 
therein). At the same time, the stellar high energy radiation fields eject electrons from heavy 
metals, and the resulting free electrons produce additional electrons by ionizing H and He atoms; 
these secondary electrons then ex cite H2 to upper levels, resulting in a chara cteristic spectrum 
of lines and continuum in the UV ([Spitzer &: Tomaskolll968l : iBergin et al.ll2004l . BO). For electron 
excitation to work efficiently, temperatures need to be high enough for neutral H to be present. The 
relatively high temperature requirements mean that the H2 detected by these means must either 
to be close to the star or to be excited by shocks. UV H2 emission has been found to be extended 
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in objects surrounded by substantial natal material, in the regions where the stellar outflow shocks 
this material, or in fast accretors, where the H2 may arise in the high density outflow itself (H06). 
However, without remnant envelopes such as the objects in this study, the only known exception 
being T Tau, the most likely place to find the required high temperatures is in the inner disk. This 
makes the UV H2 emission ideal for probing the H2 gas in the innermost regions of disks, regions 
which are difficult to access by other means. 

We obtained ACS/SBC prism spectra of a fair number of accreting Classical T Tauri stars 
(CTTS), non-accreting weak T Tauri stars (WTTS), and more evolved disks (DD), covering the 
interesting age range, ~1 - 100 Myr. Our goal was to search for UV H2 emission and study 
its evolution. The poor spectral resolution of the ACS spectra made the identification of Lya 
fluorescent lines impossible. However, we were able to identify a feature around ~ 1600 A, first 
proposed by B04 as due mostly to electron impact excitation of H2. In this letter we present and 
analyze these spectra. We show that the H2 feature is absent in all non-accreting and evolved stars 
while present in all accreting stars, and use UV fluxes to give very rough estimates of upper limits 
for the remaining surface density of H2 in the latter. 



Observations 



We obtained observations of 20 CTTS and 10 non-accreting and evolved targets using the 
Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) on the Hubble Space Telescope 
in 2007. The observations were obtained in GO programs 10810 (PI: Bergin), 10840 (PI: Calvet) 
and 11199 (PI: Hartmann). Each ACS observation consists of a brief image in the F165LP filter 
and a longer image obtained with the PR130L prism. Images appear unresolved. Offsets between 
the t arget location in the filter and prism image, including the wavelength solution, were obtained 
from luarsen ( 2006 ). The target spectrum was then extracted from a 41-pixel (1.3") wide extraction 
window. Background count rates of 0.05 - 0.1 counts s~^ were calculated from offset windows and 
subtracted from the extracted spectrum. The absolute wavelength solution was then determined 
by fitting the bright C IV A 1549 A doublet. Fluxes were calibrated from the sensitivity function 
obtained from white dwarf standard stars by iBohlinI (|2007l ) . The spectra range from 1230-1900 A 
with a 2-pixel resolution of ~ 300 at 1230 A and ~ 80 at 1600 A. 

Table [1] lists the ACS targets used in this analysis and the properties of these objects. The 
CTTS sources include 16 objects in the Taurus- Auriga molecular cloud and four sources in the 25 
Ori aggregat e in th e Orion OBla su bassociatiq n . Spe ctral types for the CTTS in Taurus are from 
Furlan et al.] i{20o d). and age s from iHartmannl (120031 ) . To correct for reddening we used the law 



towards the st ar HP 29647 ( Whittet et al. 20041 ) and estimated Ay by de-reddening the median 



photometry of lHerbst et al.l (| 19941 ) to fit the fluxes of a standard star in the region of the spectrum 
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(V to J bands) where the emission is mostly photosphericQ. We obtained accret i on lu minosities 
Lace for the Taurus sources us ing the U band exc esses following iGullbring et alj (|l998l ) , and the 
median U from photometry in iHerbst et alJ (|1994| ). The ages, sp e ctral t ypes, luminosities , yly's , 



and Lgrr. for the s ource s in 25 Ori were taken from iBriceho et alJ (|2007l ): [Hernandez et alJ (|2007l ) 



and lCalvet et al.1 ((20051) 



The non- accreting sources (WTTS/DD) were selected to have no evidence of accretion and to 
have excesses in either Spitzer Space Telescope Infrared Spectrograph (IRS) spectra or 24 and 
70 /um Multiband Imaging Photometer (MIPS) photometry, indicating the presence of debris 
disks. The sources in the TW Hydrae Ass ociation have bee n identified as WTTS by spectral 
observations which s howed Ha in emission (jWebb et al.lll999l ) and strong Li 67 07 in absorption 



(Kastner et a 
2009, . .20081 ^ ■ 



19971). The WTTS /DP and their properties were discuss e d in Carpenter et al 



Hillenbrand et al.l (120081 ). IVerrier fc EvansI (12008|), IChen et al.l ([20051) and 



Low et al 



20051 ). Examples of the ACS target spectra are shown in Figure [TJ 



We supple mented the ACS data with previously publis hed medium and high resolution STIS 



data of CTTS (Calvet et al 



2004 



Herczeg et al.ll2002l . |2004| . B04). The source properties, listed in 



Table (H were taken from I Calvet et al.l (|2004l ) for the Orion Molecular Cloud sources, and derived 
as described for t he ACS Tau r us so urces for the STIS Taurus sources . We adopt the spectral 
type and age from I Webb et al.l ( 1999 ) and Ay from 



_ ^ ■ . . , iHerczeg et al.l (120041 ) for TW Hya. Accretion 

luminosities for the STIS sample were taken from lCalvet et al.l (|2004l Umd Ingleby et al. (2009). 



3. Results 

Following B04, we identified a feature in the STIS spectra at 1600 A which is due mostly 
to electron impact H2 emission. Due to the low reso lution of the ACS spectra, we used the high 



resolution spectrum of TW Hya (jHerczeg et al.ll2004l ) to identify this feature in the ACS spectra; 



in Figure [2] we compare the feature in the observed STIS spectrum of TW Hya and in the STIS 
spectrum smoothed to the resolution of the ACS spectra. While the H2 lines are no longer observable 
in the smoothed spectrum, the feature at 1600 A is. 

In addition to electron impact H2 emission, the flux at 1600 A has contributions from accretion 
shock emission and Lya fluorescent lines (Ingleby et al. 2009). Attempting to isolate an indicator 
that is due to electron impact H2 emission, we measured the flux between 1575 and 1625 A and 
subtracted from it the continuum and the contribution from nearby strong lines (He II 1640 A and 
C IV 1550 A). Since it is unclear how strong the emission from additional sources is at 1600 A, we 
calculated the continuum in three ways. First, by joining the troughs in the spectrum on either side 
of the 1600 A feature; second, by fitting a 5th order polynomial to the entire FUV spectrum; third. 



Targets with high mass accretion rate, as DL Tau and DR Tau show significant veihng at J IjEdwards et al 
so the estimated extinction may be in error, although it is consistent with values from Taurus. 
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by adopting a continuum which assumes that the rise in the spectrum at 1600 A is due entirely 
to electron impact H2 emission. Figure [2] shows the location of the subtracted continuum for each 
method in TW Hya, and Figure [3] shows examples of the measurements for three ACS targets. 
These three methods for measuring the H2 feature luminosity were used to estimate the errors. 
Comparing the TW Hya spectra at both resolutions indicates that the feature luminosity decreases 
by ~2 in the low resolution spectrum because some of the flux is blended into the continuum. This 
error is small compared to the uncertainty in the continuum location. 

Using these procedures, we measured the luminosity of the 1600 A feature in both the ACS 
spectra and the STIS spectra smoothed to the resolution of ACS; the feature luminosities are given 
in Table 1. For the WTTS/DD, we find that the H2 feature is not observable and the values 
presented in Table 1 are upper limits based on the rms fluctuations from 1575 to 1625 A. We 
thus find that the H2 feature shows only in the accreting sources. This is not an age effect; our 
sample includes CTTS and WTTS of similar age at ~10 Myr (left panel of Figure H]) but only the 
accreting sources show the H2 feature. Moreover, we find a clear correlation of the strength of the 
feature with Lace in the CTTS (right panel of Figured]), with a Pearson correlation coefficient of 
0.68, indicating that the H2 emission d epends on the accretion properties of the source and not on 
the age. A similar result was found in ICarmona et al.l (j2007l ) , where the probability of detecting 
near-IR H2 lines was greater in sources with higher accretion rates. 



4. Discussion 



Free electrons are required for the process of electron excitation to be effective (§1). Since, in 
turn, high energy radiation fields are necessary to produce fast electrons, the absence of H2 emission 
in the WTTS/DD could in princ iple be due to a low le vel of X-ray or EUV emission in these objects 
relative to the CTTS. However. iTelleschi et al.l (|2007l ) found that there is little difference between 
the X-ray luminosities of CTTS and WTTS in their X-ray survey of pre-main sequence objects in 
Taurus. Even though there is a soft X-ray excess created in the accretion shock region of CTTS 
(Giinther et al. 2007, and references there in) , it does i iot significantly incre ase the X-ray production 
in most young stars (jTelleschi et al.ll2007l ). Similarly. iKastner et al.l (|l997l ) showed that CTTS and 
WTTS in the 10 Myr TW Hya Association have similar X-ray luminosities. Moreo ver, the X-ray 



luminosity does not decre ase significantly over the first 100 Myr of low mass stars ([Briceno et al 



19971 : iKastner et al.l 119971 ). so the CTTS and WTTS/DD in our sample should have comparable 
X-ray luminosities. 

The EUV radiation field, including emission from approximately 100 to 1000 A, is also respon- 
sible for the ionization of heavy atoms, contributing to the population of free electrons available 
to excite an H2 molecule. The EU V is difficult to investig ate because the radiation is extremely 
extincted by interstellar hydrogen. [Alexander et al.l (j2005l ) find that the EUV flux level does not 
change in the first ~10 Myr, from studies of the ratio He II 1640/CIV 1550 A. If we assume that 
the FUV level is an indicator of the strength of the EUV emission, we come to similar conclusions. 
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Figure [3] shows one CTTS and one DD that have the same FUV luminosity, so one would expect a 
strong enough EUV radiation field in both sources to create the free electrons needed to excite H2 
if it were present. However, the excess emission at 1600 A is clearly seen in the CTTS (FP Tau) 
and absent in the DD (MML 36). 

Since the high energy radiation fields in both CTTS and WTTS/DD are comparable in 
strength, the most likely explanation for the lack of H2 emission in WTTS/DD is that there is 
essentially no gas in their inner disks. Given the close relationship between the H2 feature strength 
and Lace shown in Figure [U our results suggest that H2 gas dissipates in timescales consistent with 
the cessation of accretion; when the gas is dissipated in the inner disk, there is no material left to 
accrete. 

We use the observations to make a rough estimate of the column density of H2 being collisionally 
excited. We assume that the H2 is emitted in an optically thin region of the disk with area A and 
thickness z. The emitted luminosity per unit volume is £x = hiyaxvxen'jj^, where hu is the energy 
of the emitted photon, ax the H2 cross section, v the impacting electron velocity, He the electron 
number density, Xe the electron fraction, and n^/j the number density of H2. The expected flux at 
1600 A due to electron impact excitation is then 

F16OO = ^'^^isoo^XeS^^-R^ 

where S is the surface density of H2 excited by electron impacts, mu the mass of hydrogen, R 
the radius of the emitting region, and d is the distance. In Ingleby et al. (2009) we find that the 
electron excitation model that provides the best fit to the 1600 A feature of our sample of CTTS 
with STIS spectra is characterized by a te niperature T ^ 5000 K and an electron energy of ~12 
eV . For these values, a^m n = 10~^'^cm^A ( Abgrall et al.lll997l ). According to the thermal models 
of (IMeijerink et al.ll2008l . M08), gas reaches T ~ 5000K within 1 AU of the star, which is consistent 
with the upp er limit to the extens ion of the H2 emitting region set by the STIS resolution in the case 
of TW Hya (iHerczeg et al.ll2002l ). We further assume that most electrons are capable of exciting 
H2 and adopt Xe = 5 x 10~^, as well as R ~1 AU and z ~0.1 AU (M08). Using these numbers, 
and assuming that all the flux at 1600 A is due to electron impact excitation, we get the estimates 
of S in Table 1, which for CTTS are consistent with predicted formation in the uppermost levels 
of the disk (M08). 

A similar estimate can be made for the column density of electron excited H2 in the WTTS /DD 
i n our sample, which have some dust remaining at larger radii but 110 detected IR Ho lines 
dCarpenter et al.lbood . boosi : iHillenbrand et al.lboosi : IVerrier fc Evanslboosi : Ichen et al.l[200,4 ILow et"al 
20051 ). These estimates are given in Table H We used the flux of MML 36, which is the WTTS/DD 
with the highest flux at 1600 A in our sample, to estimate the mass of H2 inside -^1 AU; we found 
that there must be less than 10~^ earth masses, 10"'^% of the MMSN, lower than the 0.01% of 



the MMSN estimated bv iPascucci et al.l (|2006l ). This has important impli cations for the forma - 
tion of terrestrial plane ts, especially if gas is needed to circularize orbits (| Agnor k. WardI |2002| ) . 
Kominami Idal (|2002l ) theorize that at least 0.01% of the MMSN must be present during the for- 
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mation of proto-planets, which form around 10 Myr according to simulations bv lKenvon &: Bromley 
(|2006l ). Our column density estimates indicate that the amount of H2 gas present in WTTS/DD 
with ages of 10-100 Myr is too small to circularize th e orbits of the terrestr ial planets being formed 
at that time. Our results support the conclusion by iPascucci et al.l (|2006l ) that there must be an 
additional source responsible for damping eccentricities, one possibility being dynamical friction 
with remaining planetesimals. Another possibility is that other species of gas exist after the H2 
has been depleted, for example, C and have been detected around the 10 Myr debris disk f3 



Pic (iFernandez et al 



2006; 



Roberge et alJ 12006 ) . C and O do not feel strong radiation pressure 



due to the low FUV flux in WTTS and therefore may remain after the H2 has been depleted 
(jRoberee et alJbood l. 
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0.3 


1.0 


1 


0.02±:g^ 




> 8.5 


CVSO 206 


K6 


0.2 


0.2 


9 


— 




> 13.9 


CVSO 35 


K7 


0.7 


0.7 


9 


o.o2±:gi 


3.2±|f 


> 16.6 


CVSO 224t 


M3 


0.1 


0.5 


9 


— 


— 


— 


OBla 1630 


M2 


1.0 


0.0 


9 






> 13.3 


STIS CTTS 
















BP Tau 


K7 


1.3 


1.0 


1 


0.23±:i 


14.1±23^ 


>41.6 


DM Tau 


Ml 


0.3 


0.6 


1 


0.08±:J? 




> 39.5 


CM Aur 


K3 


1.2 


1.1 


1 


o.i8±:?i 


l9.7±fg 


> 48.7 


LkCa 15 


K5 


1.0 


1.0 


1 


0.03±;[]i 


8.6±6| 


> 26.4 


RY Tau 


Gl 


9.6 


2.2 


1 


1 


338.0±|^[] 


> 148.4 


SU Aur 


Gl 


7.8 


0.9 


1 


o.io±:g? 




> 30.0 


T Tau 


G6 


7.8 


1.8 


1 


0.90±i6§ 


104.5±f| 


> 103.9 


CO Ori 


GO 


22.3 


2.0 


1 




303.5±ffg 


> 149.0 


EZ Ori 


G3 


5.9 


0.6 


1 


O.lOig 


20.0±f:| 


> 41.1 



Table 1 — Continued 



Object Spectral Type L Ay Age Lace H2 Feature S 

Lq mag Myr Lq lO'^L© 10"® g cm~^ 



GW Ori 


GO 


61.8 


1.3 


1 


4.7±1| 


188.2±250 


> 178.8 


P2441 


F9 


11.5 


0.4 


1 




3.4±?| 


> 29.0 


V1044 Ori 


G2 


6.7 


0.4 


1 






> 37.6 


TW Hya 


K7 


0.3 





10 


o.o3±:[jt 


2.6±3gi 


> 43.9 


ACS WTTS 
















HD 12039 


G4 







31.6 





<.021 


< 2.2 


HD 202917 


G5 


0.7 





31.6 





<.016 


< 2.8 


HD 61005 


G8 


0.6 





125.9 





<.013 


< 1.3 


HD 92945 


K2 







20 - 150 





<.006 


< 0.93 


HD 98800 


K5 


0.6 





10.0 





<.019 


< 0.84 


MML 28 


K2 




0.1 


15.8 





<.015 


< 2.4 


MML 36 


K5 




0.3 


15.8 





<.042 


< 3.3 


TWA 7 


Ml 


0.31 





10.0 





<.003 


< 1.0 


TWA 13A 


Ml 


0.18 





10.0 





<.016 


< 1.5 


TWA 13B 


Ml 


0.17 





10.0 





<.002 


< 1.5 



Note. — 

Stellar properties for binaries are from lWhite &: Ghea (j200ll ). UZ Tau A and B are themselves 
binar ies; UZ Tau A is a spectroscopic binary and UZ Tau B is a binary system (jWhite &: Ghe2 
) but is unresolved by ACS /SBC. 



2001 



CVSO 224 is a CTTS surrounded by a transitional disk ([Espaillat et al.ll2008l ) and has a very 
low M. The ACS/SBC spectrum of this target is noisy and while we cannot confidently quantify 
the H2 emission, we do see the rise in the spectrum at 1600 A which indicates its presence. 
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Fig. 1. — Sample of ACS CTTS spectra. Spectra have been corrected for reddening using the values 
of Ay listed in Tabled) Spectra have been scaled vertically for clarity. The bottom spectrum (dash- 
dotted line) in each panel is the STIS TW Hya spectrum smoothed to the resolution of the ACS 
spectra for comparison and offset by -1.2. The vertical line at 1600 A marks the center of the 
feature used to identify the H2. Left panel, from top to bottom including the offset in parenthesis: 
DP Tau (+2.7), DR Tau (+1.5), FM Tau (+0.5), FP Tau (+2.2) and GK Tau (+0.3). Middle 
panel; from top to bottom: HN Tau A (+1.9), HN Tau B (+3.0), IP Tau (+1.5), UZ Tau A (+1.2) 
and UZ Tau B (+0.45). The right panel shows ACS spectra of WTTS/DD; from top to bottom: 
HD 12039 (+3.4), HD 202917 (+2.5), HD 61005 (+2.0), HD 92945 (+1.7) and HD 98800 (+0.8). 
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Fig. 2. — Observed and convolved spectra for TW Hya. The bottom spectrum is the high resolution 
STIS FUV spectrum. The top spectrum is the TW Hya spectrum convolved to the ACS spectral 
resolution and offset by +1-0. The solid and dashed lines on the smoothed spectrum show the three 
subtracted continua. These three continua are also shown plotted on the high resolution spectrum 
and indicate that the lowest continuum may provide the best measure of the luminosity. The strong 
emission lines are labeled along with the H2 feature. 
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Fig. 3. — H2 measurements for ACS sources. The first three panels show ACS sources and the 
location of the subtracted continua, shown as the solid, dashed and dot-dashed lines. Also plotted 
are the He II and C IV emission lines as the thick solid line. The final panel compares accreting 
and non-accreting sources with the same luminosity. An excess in FP Tau is observed at 1600 A, 
which is due to electron impact H2 emission, and also between the S IV and C IV lines, which is 
likely due to blended electron impact and Lya fiuorescent lines. 
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Fig. 4. — Left: Luminosity of the H2 vs. age. Filled circles represent WTTS and open circles 
represent CTTS. For the WTTS we show only an upper limit on the luminosity of the H2. Right: 
H2 luminosity vs. Lace- The H2 luminosity is observed to increase with Lgcc- Errors on Lqcc are 
calculated using the scatter in the correlation with Ljj presented in iGullbring et al.l (|l998l ) . 



